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Theamniotic egg with its complex fetal membranes was a key innovation 
in vertebrate evolution that enabled the great diversification of reptiles, 


birds and mammals. It is debated whether these fetal membranes evolved 
in eggs on land as an adaptation to the terrestrial environment or to control 
antagonistic fetal- maternal interaction in association with extended 
embryo retention (EER). Here we report an oviparous choristodere from 
theLower Cretaceous period of northeast China. The ossification sequence 
ofthe embryo confirms that choristoderes are basal archosauromorphs. 
The discovery of oviparity in this assumed viviparous extinct clade, together 
with existing evidence, suggests that EER was the primitive reproductive 
mode in basal archosauromorphs. Phylogenetic comparative analyses on 
extant and extinct amniotes suggest that the first amniote displayed EER 
(including viviparity). 


The amniotic egg is very different from the anamniotic egg of extant 
amphibians, which lacks an eggshell and extraembryonic membranes. 
Theamniotic egg consists of a suite of fetal membranes, including the 
amnion, chorion and allantois, as well as an external shell that can be 
either strongly mineralized (as in rigid-shelled eggs) or weakly min- 
eralized (as in parchment-shelled eggs). The extraembryonic mem- 
branes enclose specific egg elements, regulate gas and fluid exchange 
between the egg and the external environment, store nutrients and 
collect waste’?. 

Where and how the fetal membranes of the amniotic egg evolved 
has been debated, and two competing hypotheses have been proposed 
(Fig. 1). The conventional, ‘terrestrial model” is that the precursor to 
amniotes laid eggs on land, similar in many respects to the directly 
developing eggs of a variety of extant amphibians, and the fetal mem- 
branes were gradually acquired so that the egg could adapt to terres- 
trial environments by retaining water inside and allowing oxygen and 
carbon dioxide to pass through the eggshell. This widely accepted 
model has been challenged by the ‘extended embryo retention model’ 
7 that the extraembryonic membranes appeared in the oviducts of 
the amniotic ancestor as specializations to control fetal-maternal 
interaction in association with extended embryo retention (EER). The 
EER model could occur with the embryo either in a post-neurula stage 


(oviparity)’ or with live bearing (viviparity)°’. Among extant amniotes, 


turtles, crocodilians and birds generally lay eggs at an early develop- 
mentalstage (non-EER oviparity), whereas most squamates (lizards and 
snakes) and mammals either display oviparity with EER or viviparity. 
Evolutionary studies based on extant amniotes give equivocal results 
about whether oviparity or viviparity arose first? ^. Circumstantial 
evidence for the EER model is the near absence of fossils of amniotic 
eggs before the Late Triassic period and the discovery of viviparity in 
many extinct amniotes as old as the Early Permian period“. Support- 
ersoftheterrestrial egg model notethat EER is absent in archelosaurs, 
including chelonians, crocodiles and birds, as well as extinct dinosaurs, 
pterosaurs and their ancestors'^"*, 

Whether the first amniote displayed EER or not is key to testing 
between the two models. As EER occurs widely among extant lizards 
and snakes (squamates) and mammals’, exploring the occurrence 
of EER among oviparous primitive archosauromorphs is decisive to 
determine the developmental stage of the first amniotic egg (Fig. 1). 
In this study, we report an articulated embryo of the choristodere 
Ikechosaurus sp. inside a parchment-shelled egg. The ossification 
sequence ofthe embryo confirms that choristoderes are basal archo- 
sauromorphs. The shell structure reveals that the aquatic choristodere 
was oviparous and presumably came ashoreto lay its eggs, like extant 
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Fig. 1 | The two competing theories for the evolution of the amniotic egg. 
Inthe terrestrial model (left), non-EER oviparity (purple) was the primitive 
condition; oviparity with EER and viviparity (blue) evolved multiple times in 
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amniotes. In the EER model (right), the evolutionarily labile reproductive mode 
of EER across oviparity to viviparity (blue) was primitive, while non-EER oviparity 
(purple) evolved multiple times in amniotes. 


sea turtles and crocodilians. The specimen, together with previous 
evidence ofviviparity inothertaxa, demonstrates that an evolutionar- 
ily labile reproductive strategy across oviparity to viviparity existed in 
choristoderes, a basal clade of archosauromorphs, and potentially also 
in other various aquatic vertebrates of the past, such as mesosaurs, 
ichthyosaurs and sauropterygians. We run phylogenetic analyses on 
extant and extinct amniotes to test whether EER and viviparity are the 
ancestral conditions in Amniota. 


Results 

Viviparity and oviparity in choristoderes 

Choristoderes are a clade of extinct diapsids that lived primarily in 
Laurasia from the Middle Jurassic period to the Early Miocene epoch 
(approximately 168-120 Myr). The gavial-like neochoristoderes were 
top predators in freshwater bodies, competing with contemporane- 
ous crocodiles”. The new specimen (MES-NJU 57004) was collected 
from yellowish white, thinly laminated tuffaceous mudstone of the 
Lower Cretaceous Jiufotang Formation (Jehol Biota, approximately 
125-120 Myr) in the Lamagou locality adjacent to Chaoyang City, 
western Liaoning, northeast China. Many choristoderes have been 
discovered in the Jiufotang Formation and the underlying Yixian Forma- 
tion in western Liaoning, including the lizard-like Monjurosuchus and 
Philydrosaurus, the long-necked Hyphalosaurus and the neochoris- 
toderan /kechosaurus. Some specimens of these choristoderes are 
associated with eggs and embryos”. 

The new specimen is a small skeleton (approximately 102.73 mm 
long; Fig. 2a) that exhibits the typical pose of a vertebrate embryo: 
curving and with the head contacting the tail". The skeleton is dor- 
soventrally flattened and exposed in ventral view covered by a thin 
layer of ferric oxide. Computed tomography (CT) scans reveal that 
the skeleton is nearly complete and all bony elements are articulated 
except for the distal end of the tail, which was slightly displaced (Fig. 3 
and Extended Data Figs. 1and 2). The embryo shows many diagnostic 
traits of /kechosaurus: the snout is long, broad and flat in front of the 


orbits and gradually tapers to about the midpoint along the snout; the 
interorbital bar is narrow; the jugal extends anteriorly to the midpoint 
of the lacrimal; the postorbital region is flared; the temporal openings 
lie largely above one another; and the parietal extends only about half 
way along the posterior edge of the upper temporal opening” (see 
the Supplementary Results for a detailed description of the skeleton). 
The fact that /kechosaurus is the only neochoristoderan in the region 
supports this assignment”. 

The phylogenetic position of Choristodera in Amniota remains 
controversial, having been placed as a basal clade of archosauro- 
morphs”, a sister group of archosauromorphs, or basal to archosau- 
romorphs and lepidosauromorphs~. The ossification sequence of the 
newembryo confirms that choristoderes are basal archosauromorphs 
(see Supplementary Results for evidence that the choristodere embryo 
isan archosauromorph). The high degree of ossification of the skeleton 
indicates that the unhatched embryo was in a late developmental 
stage. Ontogenetically, the animal was precocial or superprecocial^: 
its well-developed skull with sharp teeth suggests that it was ready to 
hunt and the relatively well-ossified pelvic girdle and hindlimbs that 
it could run and swim soon after hatching. 

We identified traces of a parchment-shelled egg around the tiny 
skeleton, which is embedded in an incomplete, oval phosphate matrix, 
demarcated by a 0.43-0.50-mm-wide halo. The halo contains mate- 
rials from both the phosphate matrix and the enclosing mudstone 
(Fig. 2d-g). The small size, embryonic pose and egg-shaped matrix 
prove that this is an embryo inside an egg. The outer edge ofthe halo 
is slightly meandering and locally folded. Loosely arranged, irregu- 
larly shaped, flake-like structures (around 50-pm thick by estimate) 
are locally present along the marginal area, surrounding pore-like 
structures (Fig. 2b,c). These features indicate that the halo is preserved 
eggshell, whichis pliable and has a thin calcareous layer composed of 
flake-like shell units and many pores”. 

There are three types of amniotic eggs: membrane-shelled; 
parchment-shelled; and rigid-shelled^?5, A mineral layer is not 
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Fig. 2| Structure of the choristodere egg (MES-NJU 57004). a, An overview of 
the embryo inside the egg. b,c, Photomicrographs show loosely arranged shell 
units and pores (black arrows) preserved along the marginal zone of the egg. 
d-g, Energy dispersive spectroscopy mapping of silica (d), calcium (e), 


manganese (f) and iron (g) showing the inferred eggshell (white triangles), which 
contains materials from both the phosphate matrix and the enclosing mudstone. 
Scale bars, 1 cm (a), 100 um (b,c), 1 mm (d-g). 


developed in membrane-shelled eggs (monotremes, a few squamates), 
isthinner than an organic layer in parchment-shelled eggs (a few che- 
lonians and most squamates) and is well developed and thicker than 
the organic layer in rigid-shelled eggs (a few squamates, most chelo- 
nians, crocodilians and birds)”. It could be argued that the described 
choristodere specimen is an incomplete viviparous egg”. In extant 
squamates, the egg of viviparous species lacks the calcified layer but is 
enveloped with a very thin organic layer (commonly less than 10 pm), 
whereas oviparous squamates lay eggs with a calcified outer layer 
and a relatively thick organic layer (usually over 30 jum)??? as in this 
specimen. A similar eggshell structure, with a thick organic layer (over 
100 um) and a very thin mineral layer (less than 10 jum), has also been 
documented in isolated eggs of the basal choristoderan Hyphalosaurus 
baitaigouensis". These shell structures reveal that the aquatic choris- 
toderes were oviparous and presumably came ashore to lay their eggs, 
like extant sea turtles and crocodilians. Putative females of the neo- 
choristoderan Champsosaurus possessed fused sacral centra and more 
robust limb bones than males, perhaps adaptations for nesting on 
land”. Other choristoderes were viviparous, such as one specimen of 
H. baitaigouensis? and Monjurosuchus splendens". The co-occurrence 
of viviparity and oviparity in one nominal species, H. baitaigouensis, 
indicates that it had a bimodal reproductive mode, that is, both vivi- 
parity and oviparity occurred in a single species as in the basalmost 
sauropsid Mesosaurus tenuidens and few extant squamates'*?. 


Macroevolutionary study 

Thenewspecimen exemplifies the complexity of parity modes in some 
early reptiles and provides information to constrain the evolution of 
reproductive strategies in archosauromorphs, reptiles and amniotes 
in general. We return to the question of whether EER is the ancestral 
condition. To test this hypothesis, we collected data on reproductive 
modes from representative taxa spanning the phylogenetic diversity of 


extant amniotes, and from all extinct taxa where information is available 
(Supplementary Table 1). The balance of data on fossil taxa inevitably 
favours those that laid rigid-shelled eggs because those are preserved 
more readily than eggs without mineralized shells. In particular, we 
could find no examples of extinct synapsids with evidence of repro- 
ductive mode. However, we compensated for this by broadly sampling 
extant taxa for which reproductive data are secure. Most oviparous 
squamates lay membrane- and parchment-shelled eggs at the limb-bud 
stage. In contrast, most amniotes that lay rigid-shelled eggs obligately 
oviposit at an early developmental stage, for instance, at the blastula 
stage in birds, the gastrula stage in chelonians and tuataras, and the 
neurula stage in crocodilians? *. In these forms, the thick calcite layer 
delays the exchange rate of respiratory gases and may prevent the 
development of the embryo before the eggs are laid^*. Therefore, both 
eggshell and developmental stage of the embryo at oviposition provide 
information constraining the ancestral state of reproductive modes. 

We coded each extant and extinct taxon for three characters: 
(1) reproduction mode: viviparous, oviparous; (2) eggshell mineraliza- 
tion: membrane-shelled, parchment, rigid; (3) EER: absent, present. EER 
was defined as amniotes that lay eggs at the limb-bud stage or later”. 
EER was identified in extinct taxa based on fossil adults that contain 
embryos atthelimb-bud stage or later or are associated with neonates, 
which were commonly identified as evidence for viviparity previously 
(Supplementary Table 1). We conservatively treated all fossils for which 
we could not judge the stage of development of an egg at oviposition 
asnon-EER. Character 2 is dependent on character 1, so viviparous taxa 
were coded as inapplicable (‘~’) for character 2. To resolve the problem 
of character dependency, characters 1and 2 were amalgamated into a 
single structured Markov model (SMM) equipped with hidden states. 
We also applied an analogous approach within a parsimony framework 
using Sankoff (cost) matrices. We conducted an exhaustive multipli- 
cative set of ancestral states analyses, accounting for: (1) different 
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Fig. 3 | CT scans of the choristodere skeleton (/kechosaurus sp.) and its 
reconstruction. a,b, Ventral (a) and dorsal (b) views of the reconstructed 
skeleton. c-f, Coronal slices showing dorsal (c) and ventral (d) sections of the 
posterior skull, proximal left forelimb and presacral axial bones (e), and pelvic 
girdle, proximal axial bones and femur (f). g,h, Transverse slices show the 
relationship of centra and neural arches in dorsal (g) and caudal (h) vertebrae. 
cdv, centra of dorsal vertebrae; ce, cervical; cl, clavicle; cv, caudal vertebra; den, 


dentary; ept, ectopterygoid; fe, femur; fi, fibula; fro, frontal; hdt, hand digit; 
hu, humerus; il, ilium; isc, ischium; ju, jugal; lac, lacrimal; mc, metacarpal; mt, 
metatarsal; mx, maxillary; na, nasal; ot, otic area; pa, parietal; pal, palatine; pdt, 
pes digit; pmx, premaxillary; po, postorbital; pof, postfrontal; prf, prefrontal; 
pt, pterygoid; pub, pubis; q, quadrate; qj, quadratojugal; r, rib; ra, radius; 

sa, surangular; sc, scapula; so, supraoccipital; sq, squamosal; st, stapes; t, 
tooth; ti, tibia; ul, ulna; v, vertebra. Scale bar, 10 mm. 


phylogenetic time-scaling methods; (2) alternative tree topologies; 
(3) exclusion of key fossils; (4) different ancestral state reconstruction 
methods, evolutionary models and optimization criteria; (5) among lin- 
eagerate heterogeneity; and (6) constraining extant node states based 
on previous studies. Together, these totalled over 100,000 individual 
analyses. Further details can be found in the Methods. 

The ancestral state reconstruction results are unequivocal for 
boththeamalgamated character (reproduction mode + eggshell min- 
eralization) and EER presence (Fig. 4). Viviparity with EER dominates 
the deeper nodes, being the most likely condition for the roots of 
Amniota (mean marginal maximum likelihood (ML)-based ancestral 
state across 100 trees for the best-fitting model and other models 
whose Akaike information criterion (AIC) difference was less than 
2 compared to the best-fitting model for viviparity-ML: 99.5-100%; 
EER-ML: 99.8-100%), Reptilia (viviparity-ML: 98.5-100%; EER-ML: 
100%), Diapsida (sensu lato) (viviparity-ML: 100%; EER-ML: 100%), 
Archelosauria (viviparity-ML: 98.5-10096; EER-ML: 99.6-100%) and 
Archosauromorpha (viviparity-ML: 99.4-100%; EER-ML: 99.8-100%), 
irrespective of which time-scaling approach or best-fitting model 
(including models whose AIC differed from the best-fitting model by 
less than 2) was considered (Supplementary Tables 3-6). 

Using a different ancestral state reconstruction also does not 
change the main conclusions: 10 of 18 maximum parsimony (MP)- 
based ancestral state reconstructions recover viviparity at the ori- 
gin of Amniota, with the remaining reconstructions mostly favour- 
ing membrane-shelled eggs as the ancestral state of amniotes (all 
MP-based reconstructions using accelerated transformation 
(ACCTRAN) favour viviparity at the origin of Amniota; Supplementary 
Tables 21-26). All MP-based ancestral state reconstructions favour 
EER as the ancestral state of Amniota. Bayesian traits (BT) results are 
consistent with ML-based ancestral state reconstructions (Supplemen- 
tary Tables 28-31 and Extended Data Fig. 4). Fixing the nodes of both 


Lepidosauria and Squamata to a non-viviparous state (Supplementary 
Tables 28-31and Extended DataFig. 5) still results in viviparity being 
the most likely condition for the origin of Amniota (mean BT-based 
ancestral state for viviparity: 87.497.796), only leading to increased 
variability in the ancestral state estimates (Reptilia viviparity-BT: 
79.8-97.2%; Diapsida (sensu lato) viviparity-BT: 72.3-99.7%) and 
potentially suggesting that viviparity re-evolved in Archosauromor- 
pha (Archelosauria viviparity-BT: 16.6-91.2%; Archosauromorpha 
viviparity-BT: 48.9-96.1%). 

Among archosaurs, the results reflect the possession of 
rigid-shelled eggs by extant birds and crocodilians, which is con- 
sistent with the proposal? that the first dinosaurian eggs were 
membrane-shelled, although the results are less clear. Most best-fitting 
models favour rigid-shelled eggs at the root of Theropoda (rigid-ML: 
81.6-98.5%; EER-ML: 0-1.5%); however, non-EER, membrane-shelled 
eggs are most likely at the root of Saurischia (membrane-shelled-ML: 
67.5-97.2%; EER-ML: 0.1-0.4%), Dinosauria (membrane-shelled-ML: 
60.5-96.7%; EER-ML: 0.1-1.4%) and probably also Archosauria 
(membrane-shelled-ML: 52-96.196; EER-ML: 0.2-13.6%). The non-EER 
results are consistent across all variant analyses but the reproduction 
mode and eggshell mineralization results are somewhat equivocal 
for later-diverging clades. The results from the best-fitting models 
are compatible with parchment-shelled eggs as the ancestral state 
for Saurischia (parchment-ML: 1-30.3%), Dinosauria (parchment-ML: 
1.1-35.6%) and Archosauria (parchment-ML: 1.3-42.5%). If we also 
consider simpler ML models with a worse fit (AIC difference greater 
than 2), where character state transition rates are more con- 
strained compared to the best-fitting models, the evidence for 
parchment-shelled eggs for these clades increases. Note, however, that 
parsimony-based analyses favour rigid-shelled eggs as the ancestral 
state of the three clades. Bayesian analyses are generally consistent 
with ML-based results but evidence for parchment-shelled eggs in 
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Fig. 4 | Phylogeny of amniotes, showing known reproduction mode and 
eggshell mineralization, and the EER of 80 extant and extinct species (tips, 
to right), and the inferred mean ancestral states for all branching points 
(larger pie charts at the nodes). The dominant inferred state towards the root 
(left) is viviparity with EER. This is a consensus tree based on a sample of 100 trees 
time-scaled using the fossilized birth-death (FBD) tip-dating method (with node 
age constraints for major clades) and component ARD model witha switch-on 
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dependency (CARD_sw) for eggshell mineralization and reproduction mode 
(best-fitting ML model based on the AIC score; Supplementary Table 6). Note that 
alternative evolutionary models do not recover viviparity as the ancestral state 

in Lepidosauria (Supplementary Tables 3-6). /kechosaurus sp. is shown in red 

and indicated by a red arrow. Elephant and bird icons reproduced from PhyloPic 
under a Creative Commons license CC BY 3.0 (elephant, T. Michael Keesey; bird, 
Chloé Schmidt). 


these three clades is generally lower; fixing the node of Lepidosauria 
and Squamata to a non-viviparous state leads to a higher variability in 
the ancestral state estimates. These discrepancies emphasize that while 
the ancestral state results are consistent for deep nodes, the results for 


later-diverging clades should be treated with more caution. A similar 
exampleis found in Lepidosauria for which parsimony andalternative 
ML and BT evolutionary models do not recover viviparity as ancestral 
(Supplementary Tables 3-26 and 28-45). 
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While the best EER model is consistently recovered as all-rates- 
different (ARD), the best model to explain character state evolution 
for reproduction mode and eggshell mineralization is more enigmatic 
and changes depending on the distribution of branch lengths, the 
time-scaling methods used and topology (Supplementary Tables 27 
and 28-45). However, the component ARD model with or without a 
switch-on dependency is always among the best-fitting ML-based 
models irrespective of which time-scaling approach or topology is con- 
sidered (Supplementary Tables 3-20). The component ARD model is 
generally also found among the best-fitting BT-based models; however, 
the evidence for alternative model hypotheses is generally lower for 
analyses with fossilized Lepidosauria and Squamata node states due 
to smaller-log Bayesian factors (BFs) (Supplementary Tables 28-45). 


Discussion 

It is widely accepted that viviparity evolved from oviparity through 
EER in squamates and mammals? and this may also have been the 
case in the various aquatic vertebrates of the past, such as choristo- 
deres, as well as mesosaurs, ichthyosaurs and some sauropterygians. 
This implies that viviparity might have derived from membrane- or 
parchment-shelled eggs with EER in these clades. Yet, our results 
strongly suggest viviparity as theancestral condition ofamniotes and 
the marine reptile clades. 

This result is subject to a number of biases: first, the phylogenetic 
placement of extinct marine reptiles within the Amniota is contested". 
Our results, however, are qualitatively unchanged if we revise the phy- 
logeny to place extinct marine reptiles within the Lepidosauromorpha 
or in other basal diapsid locations (Supplementary Tables 15-20 and 
40-45). Completely removing the extinct marine reptiles from the 
analyses has little impact on the results, with only models with a worse 
fitandsomeofthe modelsfor whichthe states of both Lepidosauria and 
Squamata were fixed a priori showing more uncertainty in the ances- 
tral state reconstructions (Supplementary Table 11-14 and 36-39). 
Further, the placement of Mesosaurus, coded ambiguously as either 
viviparous or oviparous with non-mineralized eggs” near the root of 
the Amniota does not affect these results; if itis removed, the results 
remain qualitatively the same (Supplementary Tables 7-10 and 32-35). 

Second, our dataset may be taphonomically biased. It could be 
argued that well-skeletonized embryos in adult oviducts are more 
likely to be preserved than membrane- or parchment-shelled eggs, 
andthat isolated eggs are difficult to assign to species. Consequently, 
viviparity might be overrepresented by extinct viviparous clades in 
our analysis. Given that /kechosaurus is a derived choristodere?, the 
discovery of oviparity in this presumably viviparous clade suggests 
that either oviparity evolved from viviparity or viviparity evolved 
from extinct oviparous ancestors with EER in choristoderes (the 
latter hypothesis, however, is not backed by the fossil record). This 
might also be true for extinct viviparous clades, such as mesosaurs, 
ichthyosaurs and sauropterygians that belong to basal parareptiles 
or amniotes in general and diapsids or lepidomorphs, respectively. 
As fully aquatic amniotes could not come onto land to lay their eggs 
and embryos in eggs could not survive in water, viviparity in these 
extinct amniotes must have evolved from either semi-aquatic or 
terrestrial viviparous ancestors (as implicated by our analyses), or 
oviparous ancestors that displayed EER but whose eggs have not been 
preserved in the fossil record". In either case, these early amniotes 
reflect EER as the primitive reproductive mode of amniotes (Fig. 1, 
model on the right). 

The absence of rigid-shelled eggs through the Carboniferous 
period, Permian period and most of the Triassic period? ^" has long 
been noted: if such eggs truly existed widely during this time interval, 
it could change our results entirely, pointing to rigid-shelled eggs as 
the ancestral amniote condition. However, as rigid-shelled eggs are 
more likely to be preserved than membrane- or parchment-shelled 
eggs and perhaps the partially ossified and tiny bones of retained 


embryos within the mother, this might not be a preservation bias but 
areal absence”. Conversely, given the scarcity of fossil evidence on the 
reproductive mode of Palaeozoic amniotes and the low fossilization 
potential of membrane- or parchment-shelled eggs, the possibility 
remains that the earliest amniote eggs were, indeed, membrane- or 
parchment-shelled. The addition of more extinct taxa will provide an 
important test of these results”. 

The occurrence of parchment-shelled eggs associated with EER 
in choristoderes extends the wide occurrence of this phenomenon 
among amniotes and may be the primitive reproductive mode that 
occurred before archosaurs (crocodilians, dinosaurs, birds) and 
chelonians acquired non-EER oviparity and rigid-shelled eggs'. As 
membrane- and parchment-shelled eggs associated with EER are also 
common in extant mammals and squamates, this suggests that this 
condition was common among early terrestrial amniotes, for example, 
in the first 100 Myr of their evolution from the Carboniferous to the 
Triassic, favouring a model that amniotic fetal membranes evolved 
in association with EER. Similarly, a phylogenetic study of extant 
tetrapods? inferred that many of the structures that characterize 
the Amniota (delayed deposition of eggs, large yolk mass, cellular 
yolk sac and amnion) may have evolved in an aquatic environment in 
association with delayed egg laying. 


Methods 

CT scanning and reconstruction 

The specimen was scanned at the University of Bristol, UK, on the 
Nikon X-Tek H 225 ST X-ray scanner at 225 kV and 188 uA (42.3 W) from 
a rotating tungsten target, with 2-s exposure, 1x binning, 24-dB gain 
and a3-mm copper filter, slice thickness = 48.45 um and total num- 
ber of slices = 1,142. Each scan captured 3,141 projections, with four 
frames averaged per projection. The reconstructed scan data were 
subsequently combined in VGStudio v.3 (https://www.volumegraphics. 
com). A three-dimensional (3D) model was created from the CT data 
using the segmentation tools in Avizo v.9.1.1 Lite (Visualization Science 
Group; https://www.fei.com/software/amira-avizo/). All scan data and 
3D models are available in the supplementary data. 


Scanning electron microscopy 

Samples were examined using a JEOL 8530F Hyperprobe at the School 
of Earth Sciences, University of Bristol, and a LEO 1530VP scanning 
electron microscope at the Technical Services Centre, Nanjing Insti- 
tute of Geology and Palaeontology, Chinese Academy of Sciences. 
Both instruments were equipped with a secondary electron detector, 
a back-scattered electron detector and an energy dispersive X-ray 
spectrometer. 


Phylogenetic macroevolutionary analysis 

Data collation. Data were compiled on key reproductive parameters 
for as many extinct amniotes as possible, distinguishing three char- 
acters, each with two or three states: reproduction mode: (1) vivipar- 
ity and (2) oviparity; eggshell mineralization: (1) non-mineralized, 
(2) weakly mineralized, (3) rigid; and EER: (1) absent and (2) present. 
These eggshell and parity characteristics have been documented 
widely in the literature, and we indicate exact data sources in our 
data compilation (Supplementary Table 1). We identified 59 extinct 
taxa for which eggs or viviparity had been identified (6 mammals; 
1mesosaur; 6 turtles; 13 ichthyosaurs; 5 sauropterygians; 4 squamates; 
3 choristoderes; 1 protorosaur; 2 crocodilians; 4 pterosaurs; 13 dino- 
saurs; 1 bird). We added a further 21 extant taxa, making a total of 80. 
As the basis for an initial analysis of ancestral states for Amniota and 
subclades, we compiled a supertree for the 80 taxa, using a standard 
genomic tree” as scaffold (Supplementary Table 2), supplemented by 
recent cladistic analyses of extinct groups*’. Seymouria baylorensis 
and Diadectes sideropelicus were added as outgroups to a polytomy 
including the Amniota. 
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Phylogenetic time-scaling. We obtained first and last appearance 
dates (FAD and LAD, respectively) for each taxon in our analysis using 
the Paleobiology Database*". We then time-scaled the supertree using 
four different approaches: (1) the minimum branch length (mbl) 
method”; (2) the equal branch length (equal) method**; (3) the FBD 
tip-dating method^^ * with only the root; and (4) with both the root 
and the node ages of major extant clades constrained. The mbl and 
equal methods are a posteriori time-scaling methods” that avoid 
zero-length branches by imposing a minimum branch length of1Myr 
(mbl) or by taking an equal share from preceding non-zero-length 
branches (equal). We generated 100 time-scaled trees with the time- 
PaleoPhy function of the paleotree package? with tip dates sampled 
from a uniform distribution (dateTreatment - minMax) bounded by 
FADs and LADs and the vartime argument set to 1 Myr. The FBD method 
jointly considers speciation, extinction and fossil preservation rates 
to estimate divergence times in a Bayesian framework ^ **. We applied 
aclocklesstip-dating approach using MrBayes (v.3.2.7a)**°° where the 
empty morphological matrix was generated using the createMrBayes- 
TipDatingNexus function of the paleotree package^*. Topology was 
constrained to the input supertree; uniform priors bounded by the 
respective FADs and LADs were used to calibrate the tip ages. To place 
the time-scaled trees in absolute time, M. tenuidens (Early Permian, 
Kungurian, 278.4 Myr) was used as the anchor taxon. For all node cali- 
brations, we used offset gamma distributions with a shape parameter 
of 3. Based on Benton et al.” for the root of the Amniota, we used an 
offset gamma distribution with a mean age of 325.1 Myr, a minimum 
age of 318 Myr and an s.d. of 4.099175 Myr. Additional major crown 
clade calibrations?? for the second FBD analysis were parametrized 
as follows: Mammalia (minimum age = 164.9, mean age = 182.3402, 
s.d. = 10.06911); Diapsida (255.9, 274.9603, 11.0045); Lepidosauria (238, 
245.0047, 4.044152); Squamata (168.9, 188.2463, 11.16956); Archosau- 
ria (247.1, 253.3423, 3.603972); and Aves (66, 75.91138, 5.722338). The 
default FBD and clock priors? provided by createMrBayesTipDat- 
ingNexus were kept. We disallowed sampled ancestors (prset sam- 
plestrat = fossiltip;). we ran the two FBD analyses four times, using 
four chains per run, for 1,000,000,000 generations, sampling every 
100,000th generation. We checked convergence using Tracer v.1.7.1 
(ref. 53), ascertaining an effective sample size (ESS) of all parameters 
exceeding 200 for combined traces. We used the obtainDatedPoste- 
riorTreesMrB function ofthe paleotree package to obtain a sample of 
100 time-scaled trees from the posterior, employing a burn-in of 5096. 
Before the ancestral state reconstruction, we removed the outgroup 
taxa S. baylorensis and D. sideropelicus. 


Ancestral state estimation. Hierarchical character dependencies 
have long presented a challenge in ancestral state estimation? °°. 
In our analysis, character 1 and 2 are hierarchically related; the state 
of character 2 (eggshell mineralization) is dependent on character 
1 (reproduction mode) being in a specific state (state = oviparous). 
This is equivalent to the classic tail colour problem of Maddison™. 
Recently, Tarasov? demonstrated that SMMs equipped with hidden 
states solvethe problem of modelling character complexes with hier- 
archical dependencies. In doing so, he also demonstrated the invari- 
ant nature of characters and states. These concepts are equivalent; 
characters can be transformed into states and vice versa. We followed 
the approach of Tarasov? and amalgamated characters 1 and 2 into a 
single SMM with 6 states (Extended Data Fig. 3). Ancestral states were 
estimated under two variations of this model. SMM ind assumes that 
reproduction mode and eggshell mineralization evolve independently. 
This is analogous to modelling the characters separately using two 
independent models, except for the fact that under the SMM approach, 
simultaneous changes to character 1 and character 2 are prohibited. 
Alternatively, SMM sw assumes ‘switch-on’ dependency. That is, char- 
acter 2 (eggshell mineralization) can only change state if character 1 
(reproduction mode) is in a specific state (state = oviparous). 


TheSMM ind and SMM sw approaches consider different models 
onratetransition between the character states, resulting in eight evo- 
lutionary models: a component equal-rate (ER) model (CER ind and 
CER sw:transitions between statesamong component characters share 
asinglerate parameter); a component symmetrical model (CSYM ind 
and CSYM sw: transitions between states among component char- 
acters are symmetrical); a component ARD model (CARD ind and 
CARD sw:transitions between states among component characters are 
all different); and an equal rates model (ER ind and ER sw: transitions 
between aggregated rates share a single rate parameter). The SMMs 
used in this study are summarized in Extended Data Fig. 3. 

We used an ML approach to estimate the ancestral states for the 
SMMs, applying the asr mk model function (with the optimization 
algorithm set to ‘optim’) of the castor R package”. To avoid optimiza- 
tion problems, the input trees were scaled to a tree height of 1 before 
the ML analyses. The transition matrix was fitted ten times and the 
maximum allowed number of iterations per fitting trial was set to 500. 
We used the ‘tip.priors’ argument to assign probabilities to amalga- 
mated states. M. tenuidens is a special case. Current fossil evidence 
does not enable us to determine confidently whether M. tenuidens was 
viviparous or oviparous with membrane-shelled eggs”. Again, we 
used the tip.priors argument to specify this uncertainty. For the CARD 
models, which are not time-reversible, marginal ancestral likelihoods 
were computed without rerooting the input tree. We used the AIC to 
select the best-fitting model. We then calculated the mean marginal 
ancestral states of the best-fitting model for each set of 100 input trees, 
which were plotted ona consensus tree generated using the consensus. 
edges function of the R package phytools™. Plots were generated using 
the R package strap". 

The presence and absence of EER was also modelled using asr_mk_ 
model with default parameter settings, the optimization algorithm 
set to optim, using the same input trees, and providing the two-state 
character for EER via the ‘tip_states’ argument. We used an ER model 
(EER ER) and an ARD model (EER ARD). Model selection was again car- 
ried out using the AIC. The calculation and plotting of mean marginal 
ancestral states followed the same practice as for the SMMs. 

In addition to ML, we also ran an MP-based ancestral state recon- 
struction for the same characters using the ancestral.pars function of 
the phangorn package®.. As the parsimony approach does not esti- 
mate transition rates, the evolutionary models used inthe ML approach 
cannot fully be translated into a parsimony setting. We generated 
three parsimony models for the amalgamated character (reproduc- 
tive mode and mineralization): one ACCTRAN approach that allows 
all character transitions and two most-parsimonious reconstructions 
(MPR) approaches that attempt to model either independence of 
reproduction mode and eggshell mineralization (similar to SMM ind) 
or a switch-on dependency (similar to SMM sw) using Sankoff (cost) 
matrices?^^*, The number of required steps for transitions that were for- 
bidden in the model was set to a value (step cost = 100) that would make 
it practically impossible for the transition to occur. The ACCTRAN? 
approach as implemented by ancestral.pars does not allow for cost 
matrices, thus requiring the MPR®* approach. Contrast matrices were 
used to assign probabilities to amalgamated states and account for 
uncertainty intip states. The MP-based ancestral state reconstruction 
was repeated for the EER character using ACCTRAN. 

Furthermore, we also rana BT ancestral state reconstruction using 
the package BayesTraits^^55for each set of 100 time-scaled input trees. 
As with the ML and MP approaches, the input data was formatted to 
account for uncertainty in (amalgamated) tip states. To avoid optimiza- 
tion problems, input time-scaled trees were rescaled to have a mean 
branch length of 0.01. We used the same SMM and EER used in the ML 
approach. A reverse-jump, continuous time Markov Chain Monte 
Carlo algorithm (rjMCMC) was applied to both homogeneous and 
variable rate models, the latter allowing for shifts in the rate of evolu- 
tion c? on individual branches??^?, For the models with multiple 
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transition rate parameters, we ran the multistate approach using the 
rjMCMC method with an exponential (0, 10) hyperprior. For the 
single-rate models (ER_ind, ER_sw, EERER), a uniform (0, 10) distribu- 
tion was set as the transition rate prior. Three independent MCMC 
chains per model were run for 11,000,000 iterations and parameters 
were sampled every 10,000 iterations; 1,000,000 iterations were 
discarded as burn-in. We calculated the ML of the models using the 
stepping stone sampler” implemented in BayesTraits. We sampled 
1,000 stones and used 100,000 iterations per stone. Convergence was 
assessed using the R package CODA”, ensuring that the smallest ESS 
always exceeded 200 for the combined chains. Models were compared 
using a log BF test” applied to the mean log MLs from the combined 
three MCMC chains. To calculate the log BF, the homogeneous rate 
ER sw and homogeneous rate EER ER models served as the simple 
comparison models for the SMMs and the EER models, respectively. 
Mean ancestral states were calculated across the combined three 
MCMC chains for each model and plotting followed the same practice 
used for the ML approach. Given the uncertainty in the ancestral state 
reconstructions for Lepidosauria and Squamata? "^7, we reran our 
Bayesian SMM analyses, this time fixing the nodes ofthe two clades to 
anon-viviparous state. 


Robustness tests. To test the robustness of our results, we reran all our 
analyses dropping M. tenuidens and extinct marine reptiles, respec- 
tively, from our input phylogeny. We also reran the mbl and equal 
time-scaledtrees with a modified phylogenetic position ofthe extinct 
marine reptiles, adding them either as asister taxon to Archelosauria, 
Archosauromorpha or Lepidosauria™ 8”, 


Reporting summary 
Further information on research design is available in the Nature Port- 
folio Reporting Summary linked to this article. 


Data availability 

We provide all data as supplementary data. The phylogeny used in 
this study is shown in Fig. 4. The specimen studied (MES-NJU 57004) 
is hosted at the School of Earth Sciences and Engineering, Nanjing 
University. Correspondence and requests for materials should be 
addressed to B.J. or M.J.B. 


Code availability 

All analyses in this study were conducted using readily available, pub- 
lished programs that are cited in the text. The versions of the programs 
are as follows: R v.4.1.0; ape v.5.5; castor v.1.6.7; paleotree v.3.3.25; 
phangorn v.2.7.0; phytools v.0.7-70; strap v.1.4; and BayesTraits v.4.0.0. 
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Extended Data Fig. 1| Coronal CT slices of the embryo skull. a. Approximate central view. b. Approximate middle views. Abbreviations see Fig. 1. 
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Extended Data Fig. 2| A coronal CT slice shows axial and appendicular skeletons of the embryo. Abbreviations see Fig. 1. 
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Extended Data Fig. 3 | Visualisation of the Structured Markov Models (SMMs) 
used in this study. a. Diagrammatic representation of the SMMs. SMM-ind 

= independent model; SMM-sw = switch-on dependency model. Numbers 
represent amalgamated character states. Left digit represents viviparity — 1, 
oviparity =2. Right digit represents unmineralized egg - 1, parchment egg - 2, 
rigid egg =3. b. Index matrices for the eight SMM's used in the ancestral state 
analyses. Different greyscale values represent different rate categories. Equal 


rates (ER) models: transitions between aggregated rates share a single rate 
parameter; Component equal-rates (CER) models: transitions between states 
among component characters share a single rate parameter; Component 
symmetrical (CSYM) models: transitions between states among component 
characters are symmetrical; Component all-rates-different (CARD) models: 
transitions between states among component characters are all different. 
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Extended Data Fig. 4 | Phylogeny of amniotes, showing known reproduction 
mode + eggshell mineralisation, and EER of 80 extant and extinct species 
(tips, to right), and inferred mean ancestral states for all branching points 
(larger pie charts at nodes). The dominant inferred state towards the root (left) 
is viviparity with EER. This is a consensus tree based on a sample of 100 trees 
time-scaled using the FBD method (with node age constraints for major clades) 
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and componentall rates different model without a switch-on dependency for 
eggshell mineralisation and reproduction mode (CARD ind.het) (best-fitting BT 
model based on log BF score; see Supplementary Table 31a-d), which allows for 
variable evolutionary rates on individual branches. Ikechosaurus sp. is indicated 
in red and with a red arrow. Silhouettes as in Fig. 4. 
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dependency for eggshell mineralisation and reproduction mode (CER ind.het) 
(see Supplementary Table 31d-g), which allows for variable evolutionary rates on 
individual branches. 
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